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ABSTRACT 

We present a detailed period analysis of the bright Cepheid-type variable star 
VI 154 Cygni (V=9. 1 mag, P«4.9 d) based on almost 600 days of continuous observations by 
the Kepler space telescope. The data reveal significant cycle-to-cycle fluctuations in the pul- 
sation period, indicating that classical Cepheids may not be as accurate astrophysical clocks 
as commonly believed: regardless of the specific points used to determine the O — C values, 
the cycle lengths show a scatter of 0.015-0.02 days over the 120 cycles covered by the obser- 
vations. A very slight correlation between the individual Fourier parameters and the O — C 
values was found, suggesting that the O — C variations might be due to the instability of the 
light curve shape. Random fluctuation tests revealed a linear trend up to a cycle difference 15, 
but for long term, the period remains around the mean value. We compare the measurements 
with simulated light curves that were constructed to mimic VI 154 Cyg as a perfect pulsator 
modulated only by the light travel time effect caused by low-mass companions. We show that 
the observed period jitter in VI 154 Cyg represents a serious limitation in the search for binary 
companions. While the Kepler data are accurate enough to allow the detection of planetary 
bodies in close orbits around a Cepheid, the astrophysical noise can easily hide the signal of 
the light-time effect. 

Key words: stars: variables: Cepheids - stars: individual: VI 154 Cyg - techniques: photo- 
metric - planets and satellites: detection 



1 INTRODUCTION 

Cepheids are luminous, F and G type supergiant stars, exhibiting ra- 
dial pulsations driven by the K-mechanism with periods from a few 
days up to about 100 days. Owing to their high luminosity and well- 
defined period-luminosity relations, they are primary distance indi- 
cators. The pulsation period is one of the most important param- 
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eters of a Cepheid variable and it is assumed to be stable on short 
time-scales. However, on the evolutionary time-scales, Cepheid pe- 
riods are subject to variations but these period changes become 
detectable only o ver several decades or even longer time-scales 
dSzabadosI Il983l : iTurner, Abdel-Sabour Abdel-Latif & Berdnikovl 
l200d) . 

Observations of Cepheids with ground based instruments pro- 
vide very poor light curve sampling, usually 1-2 data points per 
night which is mainly due to the long pulsation period of these 
variables. In addition, except whole Earth campaigns, observations 
usually cover a few weeks at most. The three phases of the Op- 
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tical G ra vitatio nal Lensing Exp eriment (OGLE, Soszvnski et all 
d2008allbl |201Ch ) and MACHO dAlcock et alJfT999B " projects pro- 
duced years of observations of Cepheids in the Magellanic Clouds 
but again, their light curves contain one point per night. While 
this is perfectly enough to characterize the general variability of 
the few thousand Cepheids, detailed insight into the light curve 
shape changes or short-term pe riod fluctuations is prevented by the 
daily sampling. iPoleskil J2008I) has succesfully analysed the period 
changes of Cepheids based on the OGLE and MACHO datasets and 
found variability on the scale of a few years. The near-continuous 
observations of Kepler give us a unique opportunity to study these 
variables in more detail than never before. 

Cepheids are often regarded as clockwork-prec ision, regu 



lar pu lsators an d with a few exceptions (V473 Lyr: Burki et al 
1982h. Polaris: iTurner, Savoy & Abdel-Sabour Abdel-Latif et al 



20051) : iBruntt et al.l J2008h : ISprecklev & Stevens! d2008h ) this is 



usually a valid approximation to the limits set by the sampling 
and precision of ground-based observations. State-of-the- a rt one - 
dimensio nal hydrocodes (Flor i da-Bu dapest: iKollath et alj d2002t) . 
Warsaw: ISmolec & Moskalikl {2008)) have been supporting this 
view, because after reaching their limit-cycle (in case of a single- 
mode pulsation) they are able to repeat pulsational cycles essen- 
tially forever. This of course assumes that no evolutionary or addi- 
tional processes acting on longer time-scales are included in the set 
of equations solved numerically. 

Photometric time series data obtained by space 
based telescopes h a ve b een available on a few Cepheids: 
ISprecklev & Stevensl (2008) analysed photometric data of Polaris 
observed with the S MEI instrument o n board the Coriolis space- 
craft; independently. IBruntt et al.l J2008h also st udied Polaris usin 
the SMEI and WIRE star tracker photometry; iBerdnikovl d201 
studied the periodicity of eight Cepheids based on SMEI photome- 
try; while the existing W IRE data on 5 Cephei are currently being 
analysed temndl200% . Although the precision of these data is 
superior to that of their ground based counterparts, the 0.0001 mag 
relative accuracy still hides the subtle effects to be discovered from 
the Kepler data. 

In this paper we present the first evidence that the pulsation pe- 
riod of VI 154 Cyg is jittering, i.e. it is not as stable as the present 
models predict. In Sect. 2 we briefly describe the data used, with 
details of pixel-level photometry of the saturated Kepler observa- 
tions, the methods of the analysis and simulated study about the 
timing accuracy of the Kepler data. In Sect. 3, our findings on the 
period fluctuations in VI 154 Cyg are described, while the implica- 
tions are discussed in Sect. [4] We conclude with a short summary 
in Sect. [5] 



2 DATA ANALYSIS 

The data we use in this paper were provided by the Kepler space 
telescope. The telescope was launched in March 2009 and designed 
to detect transits of Earth-sized planets. A det ailed technical de- 
scripti on of the Kepler M ission can be found in lKoch et all d2010h 



scn ption or the K epler Mls.< 
and Uenkins et al.l d2010aH bb 



The Kepler space telescope observes 105 square degree area 
of the sky in the constellations Cygnus and Lyra, and has two obser- 
vational modes, sampling data either in every 58.9 s (short-cad ence 
whose characteristics was analysed by iGilliland et all (2010)) or 
29.4 min (long-cadence), providing near-continuous time series for 
hundreds of thousands of stars. The data are divided into quarters: 
the commissioning run, Q0 (~10 d), then Ql (~34 d) followed by 



Q2, Q3, Q4, Q5, Q6, and Q7 (~90 d for each); when this paper was 
written ~600 days of observations was available to us. 

In the field of view of Keple r, there is o nly one gen- 
uine Cepheid variable so far jSzab6 etal.ll2011ab . VI 154 Cyg 
(KIC 7548061) which has been observed in long-cadence mode 
in all quarters and short-cadence mode in Ql, Q5, Q6. VI 154 Cyg 
has a mean V brightness of ~9. 1 mag and period of ~4.925 d. Pre- 
vi ous observational stu dies of this Cepheid are listed and discussed 
bv lSzabo et ai]d2011al) . 



2.1 Cepheid pixel photometry 

In lSzabo etaflfeoilah we demonstrated that the apparent large am- 
plitude variations (especially in Q2) are not intrinsic to the star, 
becau se some of the flux was lost from the assigned optimal aper- 
ture ( iBrvson et al.ll2O10l) of VI 154 Cygni. The usual output of the 
Kepler photometer is the integrated flux (light curve). However, to 
correct for this instrumental effect we have to rely on additional 
information to whic h we had no access at the time of writing the 
ISzabo etal.ld2011al) paper. 

Besides the extremely precise measurements of Kepler it is 
equally important that recently the individual pixel time-series data 
have been made availabl^ by the Kepler Team. In a number of 
cases the pixel data provided cr ucial information (1 ) on the pho- 
tocenter of the transit variations dBatalha et al]|2010l) thereby con- 
firming planet ary candidates, (2) helped to devise a custom aperture 
for RR Lyrae dSzabo et aljioich . which is too bright and is located 
very close to the edge of one of the CCDs, thereby providing a way 
to observe this importa nt prototype Blazhk o variable with Kepler, 
or (3) in a recent paper dSzab6etalj2011bl) based on pixel data we 
were able to determine which component of a close visual common 
proper motion binary is transited by a 'hot' brown dwarf or planet. 

We have downloaded all the long cadence target pixel files 
(Q0-Q7) of VI 154 Cygni. We assigned a much larger aperture than 
the original optimal aperture to sum the pixels to ensure that we 
do not loose flux, see Fig.l. It is worth noting that VI 154 Cygni is 
heavily saturated on the Kepler CCDs, so an elongated aperture is 
applied. It is also noticeable that the central (saturated) column is 
not contained in its entirety with the original optimal aperture. We 
took care to exclude additional stars captured in the downloaded 
pixels. Later it proved to be a wise decision, because by adding up 
its pixels the additional star centered at (3,13) turned out to be a 
previously unknown variable star showing frequency peaks in the 
[2,3] c/d frequency range. 

The resulting light curve is very similar to the Kepler light 
curve, except that the spurious amplitude drop has disappeared. Al- 
though light curve sections of different Qs had to be shifted verti- 
cally because of the different sensitivities of the CCDs (occurring 
due to the quarterly roll of the space telescope), we did not have 
to adjust the amplitudes which is reassuring. The only exception is 
Q0 (commissioning phase containing only two pulsational cycles), 
where inspection of the downloaded pixels confirmed that some 
flux was lost in the maximum phases, because the flux spilled out 
of the downloaded pixels, and this affected the pulsation amplitude 
and the O — C values as well. 



MAST, http://archive.stsci.edu/kepler/ 
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Figure 1. Left: original optimal aperture mask in Q2 (white), all the downloaded pixels (grey), pixels that were not downloaded (black). The axes label row 
and cloumn pixel numbers. Middle: our custom mask designed to retain all the flux. Right: Flux distribution around one of the maxima of VI 154 Cyg in Q2 
(grayscale). We deselected the pixels assigned to the near neighbour of VI 154 Cyg. 



2.2 The O - C diagram 

We studied the period sta bility of VI 15 4 Cyg applying the classical 
O — C diagram method (Sterken 2005). Thanks to the continuous 
observations of Kepler, we could determine the O — C values for 
almost all cycles, except those with gaps. For a thorough analysis, 
we used 4 different methods to calculate the O — C values that are 
described below. 

The first and second methods were to determine the times of 
maxima and minima of the light curve by fitting tenth-order poly- 
nomials around the extrema. The calculated O — C diagrams for 
the times of maxima and minima are shown in the top two panels 
of Fig.|2]plotted with black dots. The average error of the O—C val- 
ues is ±0.0021 d, which was determined by applying our method 
to a synthetic light curve (see Sect. 12. 31 ) generated using the Fourier 
parameters of VI 154 Cyg and averaging the deviations between 
the analytically derived maxima/minima and the maxima/minima 
calculated by our method. 

The third method was to measure phase shifts of each pulsa- 
tional cycle and transform this to an O — C diagram. This was done 
by fitting eighth-order Fourier-polynomial to the phase diagram of 
one pulsational cycle and used it as a master curve - allowing only 
vertical and horizontal shifts - to fit every other phase diagram. 
Combining the phase shifts and the period with the epoch, we cal- 
culated the O — C values (with an average error of ±0.0012 d which 
was adopted from the standard error of the fit of the master curve) 
which are plotted in the third top panel of Fig.|2]with black dots. 

We also determined the moments of the median brightness 
(mid value of the adjacent minimum and maximum brightness) on 
the ascending part of the light curve. Their O — C diagram are 
plotted with black dots in the bottom panel of Fig. [2] Assuming a 
0.0001 relative uncertainty in the value of both the minimum and 
maximum brightnesses of a given epoch number, the moment of 
the median brightness can be determined with an unprecedented 
precision of 0.0001 d, owing to the steepest brightness variation 
during this particular phase of the pulsation cycle. The significantly 
smaller uncertainty of the median O — C values testifies that, for 
photometric observations of low amplitude Cepheids, it is the point 
of the median brightness on the ascending branch that is the pre- 



ferred light curve feature for following the behaviour of the pulsa- 
tion period. 

For the calculation of the O — C diagrams (O stands for the 
observed, C stands for calculated times of extrema), we used the 
following ephemeris: 



C = HJD 2454969.70432 ± 4.925454 • E, 



(1) 



where C is the predicted time of the maxima (i.e. C), while E refers 
to the cycle number. The O — C values show a scatter of about 
±0.015 days (=20 min) for the full and the median and about 
±0.02 days (=30 min) for the maxima and the minima. The scatter 
of the O — C diagram is far larger than we expect from a Cepheid 
variable (see Sect.[3]and Table |2). 

We also calculated the O — C diagram using the available Ql, 
Q5 and Q6 short cadence data by applying the third method. The 
results are plotted with blue triangles in the third panel of Fig. [2] 
There is no significant difference between the resulted O — C dia- 
grams whether it is based on the long cadence and the short cadence 
data, giving a supporting evidence that the detected period jitter of 
VI 154 Cyg does not originate from the data sampling. As it is ex- 
pected, however, the average error of the short cadence O — C val- 
ues is significantly smaller than the long cadence one: ±0.0001 d. 



2.3 Timing accuracy of Kepler data 

To check the reliability of the detected O — C variations with var- 
ious methods, we designed a numerical experiment, based on syn- 
thetic data of Kepler quality, describing a hypothetic stable light 
variation (no changes in light curve or pulsation phase). The model 
was a 6th order Fourier-polynomial, fitted to the entire light curve. 
In such a way, the model described the average light curve shape of 
VI 154 Cyg. 

To produce artificial data, bootstrap of Kepler noise was added 
to the model with null signal. The noise was taken from the time 
series of VI 154 Cyg: the residuals were cycle-by-cycle fitted by 
a Fourier-polynomial, and the set of residuals were sampled with a 
replacement. Then artificial data were sent into the same algorithms 
as in the case of the observations. Since the measured distribution 
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Figure 2. The O — C diagram of VI 154 Cyg calculated with 4 different methods (black dots) using the long cadence data, while blue triangles represent the 
O — C diagram based on the short cadence data. The full description of the 4 methods is in Section |2.2| Red crosses are the resulted O — C diagrams of 
synthetic data, see in Sect. 12.31 



of O — C regards to a null-signal event, it exactly reflects the accu- 
racy of timing measurement in Kepler time series of VI 154 Cyg. 

O — C of artificial data were non-Gaussian, they exhibited sig- 
nificantly heavier wings than a normal distribution. This is because 
2-4 points are severe outliers, and they reflect probably the effect 
of uneven sampling (i.e. gaps in data). However, the method of the 
median brightness and the fit of the full light curve still results in 
standard deviations of 0.00018 and 0.00021 days, respectively, very 
small values compared to the measured effect. 

The methods of maximum and minimum brightness exhibit 
a significantly larger scatter of ~ 0.01 days, but after neglecting 
the four most distant outliers, the scatter decreases to 0.002 and 
0.003 days, respectively. This scatter is larger than in the case of 
the median brightness or fitting the full light curve, compared to 
measuring at the maximum or minimum of the light curve when the 
brightness variation is slow. However, these values are still many 
factors smaller than the detected signal. In summary, we conclude 
that the variation of O — C diagrams cannot be interpreted with 
numerical fluctuations, and the measured O — C patterns are real. 



In Fig. [3] we explain the timing method using median bright- 
ness. The upper panel shows a segment of the light curve, while 
the lower panels zoom in two regions of the same light curve. . The 
error bar of the individual data points (0.1 mmag) can still hardly 
be seen in the lower right panel, which spans 8.5 mmag in the verti- 
cal axis. Here we compare the best-fit linear function (dashed line) 
through the two consecutive data points that bracket the median 
brightness. The dotted line corresponds to the limiting case where 
the linear function is shifted vertically exactly to the lowest limits 
of the error bars, keeping the fit still consistent with the data. The 
time lag between the dashed and the dotted lines is 24 seconds; so- 
lutions with a larger time lag would misfit at least one of the data 
points, hence we take this measurement as an estimate of the accu- 
racy. 



The accuracy in timing of photometric minima and maxima is 
worse because the light curve is flat at the extrema, and the error 
bars enable somewhat larger time uncertainties. 
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Figure 3. Demonstration of the accuracy of timing the median brightness. 
The upper panel shows six cycles of the light curve. The lower left panel is a 
magnification of one individual ascending branch, and the lower right panel 
is a zoom into the region between two consecutive data points that bracket 
the median brightness. The dashed line in the small inset corresponds to the 
best-fit linear function, while the dotted line is the lowest limit that still goes 
through both errorbars. The time-lag between the two lines is an estimate 
of the accuracy, in this case 24 seconds. 



2.4 Fourier parameters 

We studied the change of the light curve shape as a possible cause 
of the scatter in the O — C diagram by examining the temporal 
variation of the Fourier parameters of the light curve. For this, we 
fitted eighth-order Fourier polynomial at the primary frequency and 
its harmonics: 



A + ^2 A i 



sin {2irift + < 



(2) 



where m is the magnitude, A is the amplitude,/is the frequency, t is 
the time of the observation, is the phase and index i runs from 1 
to 8. Then we characterised the lig ht curve shapes with the Fourier 
parameters dSimon & Teavs|[l982h . of which we show particular 
results for R21 = A2/A1, .R31 = A3/A1 as well as 02i = 02 — 
20i and 03i = 03 — 30i, 

Fig. E] shows the amplitude change and the relative variations 
of i?2i and R31, while Fig. [5] shows the relative variation of 02i 
and 03i. All of these parameters show significant cycle-to-cycle 
variations. Detailed study of these variations in the O — C and the 
Fourier parameters is described in the next section. 



3 RESULTS 

3.1 Correlations between O — C datasets 

Different determinations of O — C points show loose correlations, 
no correlation or even anticorrelat ions with each other. Results of 
a Spearman's rank correlation test dLuptonll 1993b is summarized in 
Table 1; for each compared dataset pairs, the top row shows Spear- 
man's p and the p-value is given below. This test is like a general- 
ization of the correlation coefficient: the value of p measures how 
well the dependence between two variables can be described with 
a general monotonic function. The value of p is 1 or — 1 if there 
is a perfect fit with an increasing/decreasing function, respectively; 
values in between express the relative weight of residuals which 
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Figure 5. The relative phase variations of <j>2\ (top) and cfei (bottom), with 
average errors of ±0.0013, ±0.0051, respectively. 



cannot be explained by a monotonic law. The value of p shows 
how probable is that the suspected connection between the vari- 
ables is due to a false recognition of simple numerical fluctuations. 
Id est, an intermediate value of p still safely expresses the slight 
monotonic connection between the variables if p is very small. The 
significant correlations are as follows: 

(O - C) Ma x = 0.7(2) x (O - C) Me d + 1148.6(9) 



(O - CQfuii = -0.6(1) x (O - C) Mc d + 5(5) 



(O - C*) Full = -0.52(4) x (O - C*) M ax + 600(50) 



(O - C) Max = -0.11(7) x (O - C) Mi n + 1010(90) 

(errors of the last digits are given in parentheses, absolute terms are 
in seconds). 

The most interesting is the general lack of strong dependences 
between the differently defined O — C values derived for the same 
pulsation cycles. In particular, the minimum is practically not con- 
nected to any other O — C measures; while O — C of the maxima 
shows relatively good (but still somewhat stochastic) dependence 
on the O — C of the whole cycle. Another remarkable feature is the 
anticorrelation between the O — C of the whole cycle and that of 
the median brightness or the maximum. 

This behaviour suggests two conclusions. At first, the scatter 
in O — C values is somehow reproducible, at least certain kinds of 
O — C values can be predicted to some accuracy if one knows the 
value of another O — C. Therefore, the scatter we observed in any 
O — C is not due to some unrecognized stochastic error, e.g. coming 
from the sampling and/or the evaluation of data. Perhaps the phase 
of the pulsation is changing, or there are slight local irregularities in 
the light curve. But the second finding is that the process which acts 
on the O — C values changes rapidly, this is why there are no strict 
correlations between the different determinations of O — C. These 
together suggest that the cause of the observed behaviour of O — C 
data is the irregular local fluctuations of the light curve shape. We 
shall explore this hypothesis in the next subsections. 
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Figure 4. From left to right: the amplitude (Ai) change and the relative variations of i?2i and R31 with average errors of ± 0.000042, ±0.0001, ±0.0012, 
respectively. 



Table 1. Results of a Spearman's rank correlation test of different O — C 
datasets. The sign ns means no correlation with a confidence exceeding 
85%. 
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3.2 Fourier parameters explain O — C points 

Local variations of the light curve can be described by higher-order 
Fourier-coefficients. Therefore one expects if changes in the light 
curve shape were the primary reason for O — C variations, then 
there would be connections between the Fourier parameters and the 
O — C values. If linear models are taken into account, only approx- 
imate fits are expected because the light curve shape depends on 
the Fourier parameters in a very complex and non-linear way. This 
was really observed: we generally found very slight correlations 
between the individual Fourier parameters and the O — C values. 

Another possibility is a multilinear fit of the O — C points us- 
ing the Fourier parameters. Taking the first 8 amplitude and phase 
coefficients, and additionally, the height of the minimum and the 
maximum of the individual cycles, we found that the times of the 
mean amplitude can really be predicted with 90% accuracy via a 
linear - and therefore rather heuristic - model (i.e. the scatter de- 
creases by 10% meanwhile the degree of freedom decreases by only 
17%, from 126 to 108). The accuracy was less impressive in the 
other cases, e.g. it was only 30% when fitting the times of minima 
(i.e., the residuals decreased only by 30% after the fit). 

The probable interpretation is that the O — C variations are 
really due to the instability of the light curve shape, because the 
light curve shape alone is the best predictor of O — C values. The 
existence of less impressive fits does not contradict much this inter- 
pretation, because the light curve shape is a non-linear function of 
our parameter space, and it has been foreseen that the linear fit will 
not always work. In fact, good fits can be expected only in the case 
when the parameters are combined by fortune in such way that the 
local light curve irregularities can be effectively described with a 
multilinear approximation; as it happened for the times of median 
brightness. 



Figure 6. The result of the Eddington-Plakidis test. We fit a linear trend up 
to x=15, which indicates fluctuations in the period or light curve shape. See 
text for the full description of this effect. 



3.3 Eddington-Plakidis test 

The O — C residuals (Fig.[2j were analyzed for the presence of ran- 
do m fluctuations in the pulsation period using the method described 
bv lEddington & Plakidisl (1 1 929t) 

lEddington & Plakidisl (1929) made a formalism to test for ran- 
dom, cycle-to-cycle fluctuations in period. They proposed that the 
average value (u(x)) is given as 



(u(x)) 2 = 2a 2 ± xe 2 , 



(3) 



where (u(x)} is the mean absolute difference between the O — C's 
which are x cycles apart, a is the mean observational error in mea- 
suring the time of maximum or minimum light an d e is the mean 
fluctu ation in period from one cycles to the next jPercv & Kolinl 
2000). These errors are each assumed to be accidental and uncor- 
related. Hence a plot of (u(x)) 2 against x should be a straight line 
with slope e 2 and intercept 2a 2 . 

In practice, this means that if there is a linear trend in the plot 
of (u(x)) 2 at low x (from several to several tens cycles), then it is 
assumed that random period fluctuations are present in the pulsa- 
tions of a given star and this segment of the plot is used for a linear 
fit. In other words, all wave-like fluctuations in the O — C resid- 
uals (after the removal of all trends) are assumed to be attributed 
to random fluctuations in the pulsation period or in the light curve 
shape. 



Table 2. Coefficient solutions of the Eddington-Plakidis test for four differ- 
ent methods. 



Method 


e 


e crr 


a 


Q. crr 


Full 


0.0013 


0.00054 


0.0030 


0.00085 


Med 


0.0008 


0.00042 


0.0023 


0.00066 


Min 


0.0015 


0.00058 


0.0063 


0.00141 


Max 


0.0015 


0.00048 


0.0043 


0.00105 



As it is shown in Fig. [6] a linear trend is revealed up to a cy- 
cle difference x~15. The mean random fluctuation parameters are 
summarised in Tablef2]for the four different methods. 

These results lead us to the conclusion that in short terms 
(up to a cycle difference x~15), the period and/or the light curve 
shape changes randomly. However, for long term, there is an inter- 
nal clock in VI 154 Cyg which keeps the period around the mean 
value. This feature is also supported by the fact that the O — C 
diagram of VI 154 Cygni covering 150 00 days indicates a constant 
mean pulsation period (see Fig. 13 in lSzabo et alj J201 lal) ). 

A similar Eddington-Plakidis test w as performed for SMEI 
photometric data of 8 bright Cepheids by Berdn ikovl d20irj|) . He 
found random period fluctuations, in some cases superimposed 
on period changes due to stellar evolution, but the period and its 
changes could not be followed on a time-scale of individual pulsa- 
tion cycles unlike this study of VI 154 Cyg. 



4 IMPLICATIONS OF THE RESULTS 

If the results presented here for VI 154 Cyg are common among 
Cepheids then the consequences might affect other areas of re- 
search. 

4.1 Cepheids models 

Numerical hydrodynamical models do not show such hectic pe- 
riod variations from cycle-to-cycle. This is primarily because 
these codes contain simplifi ed description of tu rbulent convec- 
tion (Florida-Budapest co de iKollath et alj j2002h . Warsaw code 
ISmolec & Moskaiikl d2008t) ~). or before the early '90s contained ra- 
diative energy transfer only. Using this type of physics the pulsa- 
tion usually reaches its limit cycle (or alternatively double-mode 
pulsation), then stays in this state essentially forever, because these 
codes are very good at conserving energy and momentum. Hence, 
the pulsation period is practically constant. Any deviation from this 
precise, regular machinery is due to choosing imperfect numerical 
scheme or coding errors. The above mentioned codes are free from 
these problems. 

Present-day models are not capab le of dealing with com- 
plex, magneto-hydrodynamic interaction jStothersl2009h that could 
cause regular or irregular modulations of the stellar structure (hence 
in the period), because of the enormous difficulties of implement- 
ing and computing 3D MHD for full stellar envelopes. However, 
Buchl er. Goupil & Kovacsl Jl993h developed a method based on the 
amplitude equation formalism to describe convective and turbulent 
random fluctuations and hence stochastic interaction between con- 
vection and pulsation. This mechanism would cause fluctuations in 
the period, but would not deform the shape of the pulsation. 

Deviations from this precise clockwor k mechanism do 
appear in certain high-luminosity models dBuchler & Kovacsl 
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1 1981 iKovacs & Buchleilll988l : iMoskalik & BuchleJ Il990l. 1 199 it 
iBuchler & Moskaiikl 19921) . namely, they can show period doubling 
or a series of period doubling bifurcations en route to chaos. These 
types of dynamical phenomena also cannot be responsible for the 
observed period variations. 

Evolutionary changes are usually happening on a much longer 
time-scale and in a regular manner, in addition the physics driving 
the evolution is usually missing from these models designed to fol- 
low only the pulsation. 

Clearly, a more precise (preferably 3-dimension) description 
of turbulent convection would be essential to model the observed 
large cycle-to-cycle period 'jitter'. 

4.2 Light curve distortion 

We can speculate that a stable global pulsation can appear as 'jit- 
tered' or distorted provided the stellar atmosphere is not quasi- 
static or spherically symmetric due to large-scale turbulent cells for 
example. This effect provides no feedback to the global pulsation, 
but perturbs the observed light curve. Changes causing local per- 
turbations of the optical depth may also induce apparent phase dis- 
turbances. This effect may be present in all cool giant variable stars 
although with different magnitude. However, the detailed modeling 
of these mechanisms are beyond the scope of the present paper. 



4.3 Searching for companions 

The usually assumed stability of the pulsation period is a fun- 
damental ingredient in finding low mass or substellar compan- 
ions around Cepheids. In light of the period jittering found in 
VI 154 Cyg, it is worth checking what kind of sensitivity can be 
expected for a perfect pulsator observed by Kepler. For this, we 
have calculated several simulations as follows. First we took a short 
subset of the long cadence light curve covering one pulsation cy- 
cle to fit a 4th order Fourier polynomial that represented the ideal 
light curve shape of the star. Then we took each time stamp of the 
combined Q0-Q7 long cadence light curve to calculate the analytic 
template, which was phase modulated by a sinusoidal term: the pe- 
riod of the modulation was set to 200 days, while the amplitude was 
set to representative values between 0.1 s and 60 s. The flux errors 
reported by the Kepler pipeline were used to add a random error 
to each point of the simulation. This way the properties of the ran- 
dom errors matched those of the original data. The simulated light 
curves were then analysed with the fourth method used in the O — C 
analysis (see Sect. 2.2), namely with the local phase determination 
using a master curve represented by a Fourier polynomial. 

In Fig.|7]we show the results for three simulations. The ampli- 
tudes of the phase modulations were 10 s, 1 s and 0.1 s, which can 
be connected to the mass ratio of the companion and the Cepheid 
using the third Kepler-law and assuming a mass for the Cepheid. 
Ignoring the sin i ambiguity of the unknown inclination angle (i.e. 
assuming an edge-on view of the system), the semi-amplitude of 
the O — C diagram depends on the Cepheid mass and the mass 
ratio as follows: 

1 f G x i/3 Ml 1/3 (l + H2) 1/3 2/3 

where Mi and Ah are the masses of the Cepheid and its compan- 
ion, respectively, P or b is the orbital period, G is the gravitational 
constant, c is the speed of light. For a 5 Mq Cepheid and a 200 days 
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Figure 7. The O — C diagram of three model light curves with sinusoidal 
phase modulations (modulation semi-amplitudes shown in the upper right 
corner in each panel). See text for details. 



orbit, the given three O — C amplitudes correspond to a companion 
mass of about 50 Mj up , 5 Mj up and 0.5 Mj up . 

The three simulated O— C diagrams in Fig.[7]clearly show that 
Kepler's precision could be enough to detect planet-sized compan- 
ions on short-period orbits around Cepheids, if the stars were per- 
fect clocks. However, the astrophysical noise in VI 154 Cyg is way 
too high to allow that kind of detections of substellar companions. 



5 SUMMARY 

The period variations for the Cepheid VI 154 Cyg have been exam- 
ined based on 600 days of photometric observations with the Ke- 
pler space telescope. In order to analyse the best quality light curve 
we performed pixel photometry on the available individual pixel 
time-series data. We calculated O — C values for all cycles using 
4 different methods: times of maxima and minima, median bright- 
ness and phase shifts of all pulsational cycles. We also determined 
the Fourier parameters of each cycle. 

We found very slight correlation between the individual 
Fourier parameters and the O — C values, suggesting that the O — C 
variations might be due to the instability of the light curve shape. 
Random fluctuation test revealed linear trend up to a cycle differ- 
ence 15 but for long term the period remains around the mean value. 

Finally, we showed that if the period jitter is common among 
other Cepheids then it needs to be taken into account in the hy- 
drodinamic models and sets a limit to detect substellar companions 
around Cepheids. 

Subtle variations in the pulsation period and shape of the light 
curve may be present in the pulsation of other Cepheids. Indi- 
rec t evidence supporting the beh aviour mentioned here was given 
by iKlaevivik & Szabados! (2009) who found that the photometric 
phase curves of some small amplitude Cepheids in our Galaxy 
show a wider scatter than expected from the measurement errors 
even if data are folded on the correct value of the pulsation period. 
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